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Abstract — In this paper, we propose novel modulation tech- 
niques using isomers as messenger molecules for nano communi- 
cation via diffusion. To evaluate achievable rate performance, we 
compare the proposed techniques with concentration-based and 
molecular-type-based methods. Analytical and numerical results 
confirm that the proposed modulation techniques achieve higher 
data transmission rate performance than conventional insulin 
based concepts. 

Index Terms — Nano communication network, molecular com- 
munication, modulation technique, isomer, diffusion, messenger 
molecule. 

I. Introduction 

Since R. Feynman's talk in 1959, people have eagerly 
pursued practical work on smaller and smaller scales, notably 
through nanotechnology (TJ. Nanotechnology, recently, pro- 
duced a new branch of research called nano communication 
networks (NCNs). NCNs interconnect several nanoscale ma- 
chines (nanomachines in short) to carry out more complex 
tasks in a cooperative manner, or to perform simple tasks (2). 
These networks are not just smaller versions of traditional 
communication networks; they have their own features and 
are applicable in many fields, including biomedical, industrial, 
military, and environmental p). 

NCNs can be realized by several methods. We can rely on 
traditional communication systems that use electromagnetic 
fields or ultrasonic wave. We have to overcome, however, 
some radio frequency (RF) device barriers (2). On the practical 
front, researchers are considering new materials, such as 
carbon-nano tubes (CNTs) and graphene |4]. Little research, 
however, has focused on channel models and human body 
absorption with terahertz bands. In addition to theses activities, 
a new concept utilizing diffusion especially for short range 
communication has been introduced (5). 

This new paradigm of communication, molecular communi- 
cation, sends/receives information-encoded molecules between 
nanomachines. One of its advantages is its biocompatibility. 
Biocompatibility is the most important and challenging issue 
confronting intra-body applications. By using biomolecules 
we can create inherently biocompatible systems that require 
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no inorganic harmful materials. Moreover, driven by chemical 
reactions, molecular communication is energy efficient. For 
these reasons, molecular communication is the main focus of 
this paper. 

The authors in (6), (7) have extensively studied the fun- 
damentals of molecular communication via diffusion. A new 
energy model has been investigated to understand how much 
energy is required to generate messenger molecules in (6) 
and concentration-based and molecular-type-based modulation 
techniques have been introduced in (7). By using a simple 
symmetric binary channel model, the achievable transmission 
rate (achievable rate hereafter) has been extensively compared. 
In analyzing their modulation techniques, however, (6), (7) 
did not clearly suggest concrete structures for the messenger 
molecules. They considered insulin-based nano networks but 
it is unclear how these could be utilized for molecular-type- 
based modulations. In this paper, to maximize the achievable 
rate with less transmit power/energy, we propose using isomers 
as messenger molecules. To consider the properties of isomers, 
we slightly modify the energy model described in (6). To 
the best of our knowledge, the proposed method is the first 
attempt to design appropriate messenger molecules for NCN 
via diffusion. 

This paper is organized as follows. Section II describes the 
channel and energy models under consideration. Section III 
explains isomers for messenger molecules and also proposes 
two modulation techniques, i) concentration-based and ii) 
molecular-type-based. We present the numerical results in 
Section IV, and Section V offers our conclusions. 



II. System Model 

For simplicity, we consider a nano communication system 
consisting of a single transmitter and a single receiver, as 
illustrated in Fig. [T] For this system, the messenger molecules 
diffuse through the medium^ (i.e. water) at body temperature. 
This paper assumes that no collisions occur among the prop- 
agating messenger molecules and/or the molecules and the 
medium. 

l ln our future work, blood will be considered as the medium to get more 
reasonable results. 



a 




Fig. 1. System model of nano communication via diffusion with a single 
transmitter and a single receiver. 



A. Channel Model 

The particles or messenger molecules released from a trans- 
mitter nanomachine spread out through the medium by Brow- 
nian motion [6]. Such motion is basically driven by diffusion, 
meaning the particles move from areas of higher concentration 
region to areas of lower concentration. A random process, 
the displacement of messenger molecules follows a normal 
distribution with zero mean. Thus, the standard deviation a of 
the displacement can be obtained as follows: 
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where, c indicates the concentration of Brownian particles at 
time t at point x, and D represents the diffusion coefficient 
of the particles calculated from the Boltzmann constant CK&), 
temperature (T), the viscosity of the medium (rj) and the 
radius of a messenger molecule (r mm ). Eq. ([I]) is Fick's 
second law of diffusion. By solving this partial differential 
equation, we obtain the general equation for the concentration 
of the particles represented in ([2]). The first moment of ^ 
is zero, which indicates the displacement has a zero mean, 
and the second moment becomes the variance. Hence, the 
displacement has a standard deviation of \f2Dt, and finally 
we have 



AX ~AT(0,2Dt). 



(4) 



When n messenger molecules are transmitted by the trans- 
mitter nanomachine, the molecules have a probability of hit- 
ting the receiver nanomachine. We represent this as a binomial 
distribution of n times of trials with a probability of Phu for 
each. Phu is determined by the symbol duration (T s ) and the 
distance between the transmitter and the receiver (d), which are 
both affected by the diffusion coefficient. If n is large enough, 
and nPhu is not zero, binomial can be approximated as a 
normal distribution. In addition, we have to take into account 
the overflow molecules from the previous symbol (6). The 



total number of molecules during a symbol duration becomes 
different depending on the current and previous symbol. 



Binomial(n, P hit (d, T s )), 
N c ~ Af(npi,npi(l - pi)), 
N p ~ Nnp 2 ,np 2 (l -P2)) ~ 



N(npi,npi(l - pi)) 



where, N c denotes the number of molecules transmitted and 
received during the current symbol duration, and N p is the 
number of molecules transmitted from the previous symbol 
duration but received during the current symbol duration. For 
simplicity, we use p\ for P^u during T s and p 2 for Put during 
2T S . 



B. Energy Model 

We use an energy model to that described in [6]. To prevent 
them from interacting with others during propagation, the 
messenger molecules are encapsulated in vesicles. The vesicles 
are then carried to the boundary of the machine (e.g. the 
eukaryotic cell), and released into the propagation medium. 
Here we use r unit to denote the machine radius size. 

The energy needed in these steps are calculated. Moreover, 
we can obtain the cost of synthesizing messenger molecules 
(hexoses as an example) by the enthalpy of formation AH. 
Then the final energy model can then be presented as fol- 
lows (6): 
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where, E T is the total energy cost required to transmit n num- 
ber of molecules, Es is the synthesizing cost of one hexose 
molecule calculated from the sum of bond energy (e.g., the 
enthalpy change), Ey is the vesicle- synthesizing cost having 
a radius of r v . Ec is the cost of intra-cellular transportation, 
Ee is for membrane fusion, and c v is the capacity of one 
vesicle that is related to the radius of messenger molecules 



III. Modulation Techniques 

In molecular communication via diffusion, the unique prop- 
erties of the messenger molecules can determine modulation 
techniques. Two such techniques that have been proposed 
include the use of concentration and type of messenger 
molecules |7|. A usable messenger molecule itself, however, 
has yet to be proposed. Thus this paper proposes practical 
messenger molecules, as well as analyzes and compares their 
achievable rates by applying different modulation techniques. 
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Fig. 2. D-form isomers of aldohexoses. 
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Fig. 3. Mirror image of D- and L- glucose. They are enantiomers. 

A. Isomers for Messenger Molecules 

The most important thing we have to be careful when 
designing messenger molecules is that has to be non-toxic 
to human body. The messenger molecule suggested in (6j, 
however, is highly flammable (i.e., hydrofluorocarbon), which 
means it may be inappropriate as a messenger molecule. 

For several reasons, potential candidates for messenger 
molecules are isomers, molecules with the same number and 
types of atoms [8]. First of all, they consist of the same type 
of atoms, lightening the burden on the transmitter nanoma- 
chine synthesizing the messengers. For numerical analysis, 
this paper uses the isomers known as hexoses, especially 
aldohexoses. Note that we can select from the aldoses family, 
hexoses, pentoses, tetroses, or trioses based on the required 
modulation order. 

Aldohexoses are monosaccharides with the chemical for- 
mula CqH^Oq. They have four chiral carbon atoms that give 
them 16 (=2 4 ) stereoisomers. Fig. [2] illustrates eight kinds of 
D- form diastereomers. The enantiomers of each molecule 
are another set of eight L- form diastereomers. Therefore, 
aldohexoses have 16 different shapes. Here, diastereomers are 
stereoisomers that are not enantiomers (mirror-image isomers). 
For example, D-glucose and L-glucose as shown in Fig. [3] 
are enantiomers. D-glucose and D-galactose, isomers but not 
mirror images, are diastereomers. 

When each isomer is dissolved in aqueous solution, it 
mostly exists as a cyclic form. D-glucose, for instance, 
undergoes nucleophilic addition reaction generating four 
cyclic anomers, a-, f3- forms of D- glucopyranose and D- 
glucofuranose (9). We consider, however, only two pyranose 
forms since they predominate with 36 and 64 percentages. If 
the functional groups attached to carbon number 1 (CI) and 
the C5 shown in Fig. [4] have a trans -structure, it is called 



a- form, and if a cis-structure, (3- form. They interconvert 
each other in solution through a process called mutarotation. 
Therefore, by deploying hexoses group in the system, there 
are a total of 32 different isomers. This means that we can 
increase a modulation order up to 32, i.e., 5 bits per symbol. 

B. Molecular Concentration Based 

When the concentration of messenger molecules is used, 
the technique is known as concentration-based modulation, 
originally introduced in (7). It is also referred to as CSK 
(concentration shift keying). In binary CSK (BCSK), one 
threshold exists, and a receiver nanomachine decodes the 
symbol as ' 1 ' if the number of received messenger molecules 
exceeds the threshold, '0' otherwise (6). Generally, 2 n -CSK 
systems transmit n bits per symbol, and requires 2 n -l number 
of thresholds. In this system, there is, theoretically, no limit 
in the modulation order. As the modulation order increases, 
however, so does the error probability since the minimum 
distance between two neighboring thresholds decreases. It only 
uses one kind of molecule, and D-glucopyranose is used for 
analysis here. Below is a probabilistic analysis for BCSK 
systems: 

p a (o,o) = n(o,o,o) + n(i,o,o) 

= \iP(N n <T) + P(N p + N n <r)} 
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where, P a (X,Y) indicates the probability of X sent and Y 
received, and Pb(Z,X,Y) is the probability of X sent, Y 
received, and Z previously sent. Q(-) is the tail probability of 
the normal distribution. 

C. Molecular Type Based 

When different types of molecules represent different sym- 
bols, the technique is known as molecular- type-based modu- 
lation, referred to as molecule shift keying (MoSK). Unlike 
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TABLE I 
Simulation Parameters. 





Fig. 4. D-glucose and its a- and (3- anomers. They undergo rapid 
interconversion in water, which is called mutarotation. 



the work in (7), we here use a set of isomers. We thus 
name it IMoSK (isomer-based MoSK). IMoSK requires only 
one threshold, making it simpler than the CSK system. For 
systematic analysis, we can choose one among the several 
aforementioned isomer sets; a modulation order can be de- 
termined by the sets used. For example, hexoses have a 
modulation order of up to 32, trioses of 4. For simplicity, 
we apply, as did in (7), the addictive white Gaussian noise 
(AWGN) model. The mutarotation effect is also considered. 

1) B-IMoSK-AWGN: Only AWGN considered, 



P a (a, a) awgn = Pb{<x, a, a) + P&(/3, a, a) 

= \[P(N p + N c + N n > r)P(N n < r) 
+ P(N C + N n > t)P(N p + N n < r)], 

P a (ot, P)awgn = Pb(ot, a, /3) + P b (fi, a, 0) 
1 
4 

+ P(7V p + 7V n > r)P(7V c + 7V n < r)]. 



[P(7V n > r )P(7V p + N c + 7V n < r ) 



Note that P a (/?,/?) and P a (/3,a) can also be calculated 
similarly. Due to the page limit, we omit it here. 

2) B-IMoSK-Mutarotation Considered (a-D-glucopyranose 
<^> f3-D-glucopyranose): When a- and /3- D-glucopyranose 
are chosen, there is a possibility of incorrect decoding, such 
as a- sent, f3- received, or f3- sent, a- received due to 
mutarotation. Thus we derive an error probability considering 
the mutarotation process. The a- or j3- form sent varies its 
number with time fTO) , and the number can be calculated by 
observing the specific optical rotation. R ta or R t p (observed 
specific optical rotation at time t) minus R eq (at equilibrium) 
divided by R a or Rp (at time 0) minus R eq has a linear 
relationship with time (T s ) at about 36.5° C. R ta and P t /3 
values are calculated by R a and Rp, and R eq cab be found 
in (TO). The number of a- and (j form existed after time T s 
is obtained as below. If the number of /?- form exceeds the 
threshold value after T s when a- sent, - value is added to the 



Parameters 


Values 


Phit for T s 


0.6097 


P hlt for 2T S 


0.7208 


Radius of the hexoses 1 12] 


0.38 nm 


D of hexoses J5J 


597.25 iirn 2 /sec 


^-f^hexose 


1211 kJ/mol 


T s 16] 


5.9sec 


Viscosity of the water 


O.OOlkg/sec • m 


Temperature (body temperature) 


36.5 °C = 310K 



error term. Thus, P a (a,f3) can be calculated as 



Rt, 



Rp 



Rta 



R, 



eq 

aR n 



0.99 
~3600 l 

■PR? 



aR a + (n — a)Rp 



(Rt 



a + P 
- Rp)n 



If /3 > r, P a (a, /3) = P a (a, /^Wgn 
where, P P „ = 52.7°, Pv=112.2°, and P 



n 

18.7° ifTOl 



The 



parameters a and /3 indicate the number of a- form and f3- 
form molecules, respectively, and n is the number of total 
molecules transmitted. Optical specific rotation of the chemical 
compound is defined as the observed angle of optical rotation 
when plane-polarized light passes through PJ-glucopyranoses^] 

3) 32-IMoSK: When using hexoses, the system has the 
maximum modulation order of 32. Probabilities of X sent and 
Y received for all X and Y values are obtained similarly and 
thus omitted here (see p3|). 

IV. Numerical Results 

Assume that all hexoses have the same physical properties; 
size, diffusion coefficient, and enthalpy of formation. P^u 
value is calculated by the same numerical calculation used 
in (6). We approximate the value for the hexoses by assuming 
it to be proportional to the diffusion coefficient. Therefore, 
we have the hitting probabilities shown in Table |T| and apply 
these into the proposed scheme explained in Section [Hi] Here, 
we define the achievable rate R that maximizes the mutual 
information I(X;Y) as follows 

P(X,Y) 
P(X)P(Y)' (5 ) 

P = max/(X ; y) 



I(X;Y) = Y P(X,Y) \og 2 



where, P(X,Y) denotes the joint probability of X and Y 
while P(X) and P(Y) are the probabilities of events X and 
Y. In d5|, the threshold values vary from 1 to 1000. 

Fig. [5] compares the achievable rate of the proposed method 
using hexoses with the conventional insulin based method in 
BCSK and B-IMoSK systems. It is remarkable that we obtain 
about 8 dB of SNR gain in both BCSK and B-IMoSK systems. 

2 It could be measured by a polarimeter, and there is a linear relationship 
between the observed rotation and the concentration of the compound [H] 
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Fig. 5. Achievable rate comparisons of the conventional insulin based and 
the proposed methods using two kinds of aldohexose isomers. 



Fig. 7. Achievable rate comparisons. Triose and hexose used. 
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Fig. 6. Achievable rates of 32-IMoSK and B-IMoSK. Thirty two kinds of 
isomers of aldohexose used for 32-IMoSK. 



In addition, the B-IMoSK system shows a better data rate 
performance than BCSK. This is due mainly to the size of the 
proposed messenger molecule being much smaller than that 
of insulin. Hence, the transmit energy of the proposed method 
is much less than that of the conventional method. This result 
will be discussed in more detail in (T3j. 

Fig. [6] shows the achievable rates of 32-IMoSK and B- 
IMoSK using hexoses as messenger molecules. As can be 
seen from Fig. [6] the 32-IMoSK system has the maximum 
rate of 5 (bits per symbol), and the B-IMoSK of 1 (bit per 
symbol). In Fig. [7] we compare the achievable rates using 
trioses and hexoses. Obviously, trioses achieve higher SNR 
gain than hexoses due to their smaller sizes. Trioses, however, 
have the data transmission limit 1 bit per symbol. From this 
result, we can conclude that trioses can be selected for a low 
data rate system with a higher transmission reliability and 
hexoses for a high data rate system. 

V. Conclusion 

To make nano communication feasible in practice, this 
work proposed novel modulation techniques using isomers as 



messenger molecules. We first introduced energy and channel 
models for our system. Next we proposed several modulation 
methods able to support up to 5 bits per symbol. We also com- 
pared the achievable rate performance with existing modula- 
tion methods (concentration-based and molecular-type-based). 
This work differs from prior work in that it proposes practical 
messenger molecules and provides guidelines for selecting 
from among several possible candidates. Future work, we will 
consider, to achieve more modulation degrees of freedom, 
the ratio of enatiomers. To increase the transmission data 
rate further, we will also consider multiple sets of messenger 
molecules. 
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